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We present a calculation for the polarizations of J/tp and tp' produced in 
the hadron collisions at the Fermilab Tevatron. Various color octet channels 
including ^Sq , ^Pj , and ^S[ as well as contributions from Xcj decays are 

X" 

Vh ' considered in the kt factorization approach. We find that in a rather wide 

a ■ 

range of the transverse momenta of J ftp and tp' , the production rates could 

be dominated by the ^^q channel, and the predicted polarizations from the 

^5*0 channel and Xcj feeddown contributions are roughly compatible with the 

preliminary CDF data. This might provide a possible release from the conflict 

between the NRQCD collinear parton model calculations and the CDF data. 
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Heavy quarkonium production in high energy colhsions provide important information on 
both perturbative and nonperturbative QCD. In recent years, heavy quarkonium production 
has attracted much attention from both theory and experiment. To explain the J /ip and ip' 
surplus problem of large transverse momentum production at the Fermilab Tevatron |I| , the 
color-octet production mechanism was introduced for the description of heavy quarkonium 
production |^ based on the NRQCD factorization framework |^. However, most recently 
the CDF collaboration has reported preliminary measurements on the polarizations of the 
promptly produced charmonium states [Q, which appear not to support the color-octet 
predictions that the directly produced S-wave quarkonia have transverse polarizations at 
large px HiH- In 0> ^^^ authors further considered the feeddown contributions from Xc 
decays, and found the prompt J /ip polarization disagree with the CDF data at large pr by 
3 standard deviations. In p[, we have performed a calculation of the production rates for 
J/ip and ip' at Tevatron in the newly advocated kt factorization approach 0, and we find 
that within the color-singlet model the production rates of J /ip and ip' can be enhanced 
by a factor of 20 compared to the leading order predictions in the naive collinear parton 
model, but are still below the data by at least one order of magnitude. Therefore in order 
to explain the J /ip and ip' productions at Tevatron we still need to call for the contributions 
from color-octet channels. This conclusion has been confirmed by another independent study 
on the J /ip production in the kt factorization approach [0. In this letter we will report a 
calculation for the contributions of color-octet channels not only to the cross sections but 
also to the polarizations of the J/ip and ip' . 

The /c(-factorization approach differs greatly from the conventional collinear approxima- 
tion because it takes the non-vanishing transverse momenta of the scattering partons into 
account. This approach was previously advocated for the study of heavy quark produc- 
tion at high energy hadron colliders |TT| , p!2| . Recently, more phenomenological studies were 
performed along this direction by several groups [p!3Hl5|, and some successes seem to be 
achieved by comparing the calculations with experimental data on heavy quark production 
at the DESY HERA and the Fermilab Tevatron [[l^ -[T6[|. Especially for h and h production. 



the kt factorization approach may provides a successful resolution for the long standing con- 
flict between the NLO collinear parton model prediction and the experimental data at the 



Tevatron |T5[. From these and the most recent studies on charmonium production 
we see that for heavy quark and quarkonium production in the present pt region at the 
Tevatron, the improved kf factorization approach may play a very important role and bring 
significant enhancements compared with the conventional collinear parton model. For the 
b quark production the enhancement is about a factor of 2, while for the color-singlet J/ijj 
production a more dramatical enhancement, which is about a factor of 20, has been found 
P]. In this paper we will work along this line, to further study J ftp polarization at the Teva- 
tron by including the color-octet channels. We find that the kt factorization approach may 
provide a possible way to understand both the observed production rates and polarizations 
of J/ip and ip'. 

Under the kt factorization approach, a process can be factorized as the hard scattering 
cross section convoluted with the so called unintegrated gluon distribution of the proton. 
For example, the polarized cross section for a charmonium state H at the Tevatron can be 
written as the following form, 

— — = / dxidx2d qird gaT o 2 ^ ^{9i92 ^ H^ ' + X), (1) 

where Pt and y are the transverse momentum and rapidity of the charmonium state H, and 
A denotes its helicity. /(x; q^) is the unintegrated gluon distribution function, which can be 
related to the conventional gluon distribution by 



M 



dq^ 



xg{x,fi)= -^f{x]qj). (2) 

J qj, 

In Eq. (1), the hard scattering cross section cr{H^'^^) describes the production of H from 
two off-shell gluons fusion process, which will depend on the incident gluons' longitudinal 
momentum fractions Xi, X2 and transverse momenta gir, q2T- In the collinear approxima- 
tion, the cross section a does not depend on qiT and q2T, and then Eq. (1) will go back to 
the conventional collinear parton model result for the production of state H. However, as 



mentioned above, the incident partons' qt may play an important role for heavy quarkonium 
production, and may not be neglected. To derive the cross section o", we closely follow the 
BFKL approach in which gauge invariance of the interaction vertices is guaranteed order 
by order. For example, to get the gauge invariant amplitude for J /if) production from ^5*} 
channel, we must include extra term contributions (see Fig. 1(c) of |Q]) which will finally 
result in an effective (Lipatov) vertex for three-gluon interaction. 

For direct J /ip production at the Tevatron, because the color-singlet contribution under- 
estimates the experimental data p| by at least an order of magnitude in the kt factorization 
approach, we need to consider the color-octet contributions, for which we will have the LO 
(leading order) gluon fusion processes, 

gg^cc['Sf\'pf\'Sf^^J/i, + X, (3) 

where various possible color-octet channels for the cc are considered. Note that these 2 — i> 1 
subprocesses do not contribute to J /ip {ip') production at large px in the collinear parton 
model. However, in the kt factorization approach, these processes contribute to large pt 
production, because the initial partons' transverse momenta will lead to J/ip having trans- 
verse momentum. In this context, the transverse momentum of J/ip will be the sum of the 
transverse momenta of the two incident partons, -Pf = Qit + q2T- For Xcj production, the 
inclusive production rates have been calculated in [^], where they find there is no need for 
the color-octet contributions to Xcj production. So, in this paper we will also only consider 
the color-singlet contributions to Xcj polarizations, and then calculate their feeddown con- 
tributions to J/ip polarization. For numerical calculations, we set rric = 1.5GeV and choose 
the unintegrated gluon distribution of [|1^]. We set the scales /i^ for the strong coupling 
constant ^^(/i^) in the hard scattering cross section a to be qfrp for the interaction vertex 
associated with the incident gluon gi, and g|j. for the vertex associated with q2 [p|,p|,p!8|. 



We first give the numerical results for the Pt distributions of the inclusive cross sec- 
tions for direct J/ip and ip' production at the Tevatron compared with the experimen- 
tal data. Fig. 1 is for J/ip, and Fig. 2 for ip'. As a qualitative estimate we will use 



just one individual color-octet channel to fit the cross sections observed by CDF. The 
solid lines correspond to the ^5*0 channel by taking the color-octet matrix elements as: 
{Oti^So)) = O.lOGel^^ fo^ j/^ ^^j (of (i^o)) = OMGeV'^ for f, while other matrix ele- 
ments are set to be zero, {Og (^-Po)) = (Cg ('^'S'l)) = 0. The dashed lines correspond to 
the 3pjs) channel: {Oti^Po))/ml = 0.025GeV^ for J/i/j and {Ot'CPo))/ml = 0.005GeV^ 
for ^', while {Of^'\^So)) = {Of'^'\^Si)) = 0. The dotted lines correspond to the ^S?^ 
channel: (C^(^^i)) = O.OOTGel/^ for J/ijj and {OfCSi)) = 0m25GeV^ for ^', while 
{Ot^^'\^So)) = {Ot^'^'\^Po)) = 0. For J/iP production, from Fig. 3 we can see that the 
^5*0 and '^PJ channels can both rather well describe the Tevatron data, especially for the 
shapes of the differential cross sections. Whereas the ^S"} channel gives an inadequate shape 
for the differential cross sections compared with the experimental data. For ip' production, 
at first sight all the three curves seem to roughly agree with data. However, we can still 
see that the ^S*} channel overestimates the production rate of ip' at large Pt- For both 
J/-ip and -ip' the theoretical predictions for ^5*} channel drop more slowly as Pt increases 
than the experimental data. This is mainly because for ^S"} channel there are additional 
terms (see Fig. lb and c of [^), which give main contributions at large Pt, but these terms 
do not contribute to the ^Sq and ^PJ processes. The disagreements between theoretical 
predictions for the ^S"} contributions and experimental data show that in the present pt 
region {pt < 20GeV) at the Tevatron, the ^5*1 channel may not be the dominant one to 
the direct J/ip and ip' production in the kt factorization approach. This result is very dif- 
ferent from the collinear parton model calculations, where the ^S*} channel overwhelmingly 
dominates the S'-wave quarkonium production at moderate and large transverse momentum 
at the Tevatron P,|^. 

This difference has significant consequence on the polarization predictions of J/ip and 
ijj'. We recall that in the collinear parton model the NRQCD formalism predicts S'-wave 
quarkonia being transversely polarized dominantly due to the ^5*} gluon fragmentation 
contributions |§J§]. However, in the kt factorization approach, the ^S*} channel may not 
dominate the S'-wave quarkonium production, so the previous conclusion of transverse po- 



larization at large Pt may be changed. In the following we will study the polarizations of 
J/ip and ip' in the kt factorization approach. In Fig. 3, we present our result for the tp' 
polarization at Tevatron. Because it receives no feeddown contributions from higher lying 
charmonium states, the polarization of ip' is only due to the direct production from the 
color-octet processes. In this figure, we give each individual contributions from these octet 
processes in the kt factorization approach. The solid line is for the ^Sq contribution, the 
dashed line for the ^PJ , and the dotted-dashed line for the ^S"} . From this figure we can 
see that the ^S*} channel contributes to the transversely polarized ip' as it contributes in the 
collinear parton model approach. Whereas the ^5*0 channel contributes to the unpolarized, 
and ^Pj contributes to the longitudinally polarized. If ^S'q could dominate ip' production 
at the Tevatron (see Fig. 2), the ip' will be unpolarized. This is consistent with the exper- 
imental data, which show that in the available pt region the ip' seems to be unpolarized 
(though with large errors). So, the experimental data on ip' polarization support the ^Sq 
dominance on tp' production. 

We then study the J /ip polarization. The polarization data of J /ip are for the prompt 
production, i.e., including the direct production and the feeddown contributions from Xcj 
and ip' decays. The ip' feeddown contributions to J /ip polarization can be obtained from 
the above ip' polarization calculations multiplied by the decay branching ratio of ip' to 
J/ip. (Note that the observed transitions of ip' to J/ip involve no spin flips, so this part of 
contribution has the same behavior as ip' polarization.) To obtain the feeddown contributions 
from the Xcj decays, we must convolute the Xcj polarized cross sections with their decay 
branching ratios of Xcj -^ J l^^l- We find that in the kt factorization approach XcJ feeddown 
contributes to J jip being a little transversely polarized. The polarization parameter a for 
J jip from XcJ decays approaches to 0.5 at large p^. This is different from the collinear parton 
model calculations which predict that a approaches to about 0.24 at large p^ 0)0 • (Note 



that the polarizations of XcJ themselves are very interesting for testing the kt factorization 
approach, and has been presented in ||2^). All of these contributions including the direct 
contributions and the feeddown contributions to the J jip polarization are plotted in Fig. 4. 
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As in Fig. 3 for ip' , we also plot in this figure the individual contributions from the three 
color-octet processes for the direct polarizations, and then combine them with the Xcj and 
ip' feeddown contributions for the prompt polarizations compared with the experimental 
data from the CDF collaboration. The thinner lines are for the direct polarizations, and 
the thicker lines for the prompt polarizations. The solid lines are for ^5*0 contributions, 
the dashed lines for ^PJ , the dotted-dashed lines for ^S*} , and the dotted line for the Xcj 
feeddown contribution. We see that in general the polarizations increase as pt increases. 
For example, for the ^Pj direct production the value of polarization parameter a{ilj) (being 
longitudinal) varies from —0.63 at pr = 5GeV to —0.38 at pr = 20GeV. For the ^S*} direct 
production q;(V') (being transverse) varies from 0.72 to 0.96 in the same Pt range. The Xc 
feed down contribution will reduce the extremities of a{ilj) in direct productions. From this 
figure we find again that the ^5*0 dominance together with the Xc feed down contribution 
will lead to the predicted prompt J/ip polarization in rough agreement with the Tevatron 
data except the one at about px = lQ.7GeV which, however, has very large errors. We 
note that this seemingly nontrivial behavior of J/ip polarization at px = 16.7GeV (being 
longitudinal) can not be explained consistently with the moderate Pt data (being transverse) 
in our approach. If the CDF data at pr = 16.7GeV is further confirmed by CDF or DO 
experiments in the future with higher statistics, the above mechanism may fail. 

Concluding our analysis, we find that the ^5*0 channel dominance is the only reasonable 
realization of the polarization data of J/ip and ip' in the /c^ factorization approach. To see 
this point more clearly, we can perform a x^ analysis for theoretical comparison with the 
experimental data. For ^' polarization data, ^5*0 channel gives the smallest value of x^ 
{x'^/dof = 0.61) as expected. Moreover, inclusion of other channels does not improve the 
fit. For example, including the ^S*} contribution would make the fit even worse, and when 
including ^PJ contribution the fit is in the same manner as ^5*0 alone but with a little 
higher x^ /dof (= 0.83). For J /ip polarization, the result for the fit is similar. I.e., ^5*0 
alone gives the smallest x^ value {x'^/dof = 1.85), and inclusion of ^PJ gives a even worse 
fit, while inclusion of ^S*} gives the same quality fit as ^5*0 alone with a little higher x^ /dof 

7 



(= 2.13). So, with the available statistics at present, the polarization data of J/ip and ip' 
may indicate a strong signal of ^5*0 dominance for J /if) and if)' production. 

The ^5*0 dominance has a straightforward consequence that the other two matrix ele- 
ments (pfi^Si)) and (Cg (^Po)) must be much smaller than {Ofi^So)). This is quite different 
from the naive expectation based on the NRQCD velocity scaling rules, which predict no 
big difference between these three matrix elements ((Og (^5*0)) scales as f^, while {Ofi^Si)) 
and (Og (^Po)) both scale as v^). But, the values of (Cg ("'^'S'o)) we used in Figs. 1-2 are 
smaller than the relevant color-singlet matrix elements by about an order of magnitude (e.g., 
the observed J /if) leptonic decay width leading to {Oii^Si)) = l.^QGeV^). This is roughly 
consistent with the NRQCD velocity scaling rule estimate. As for the other two matrix el- 
ements (Og (•^S'l)) and (Cg (^Pq)), there might be some dynamical reasons (associated with 
the nonperturbative evolution processes) to suppress their values and violate the velocity 
scaling rules pl|, or there might be some new counting rules for charmonium system to 



suppress their values [22 



Some theoretical uncertainties must be addressed before concluding. First there is an 
uncertainty coming from the scale choice used for the numerical calculations. We have set 
the scales to be q^j, or g^r which are also adopted by some previous studies [|I^,|rH] . If we set 



these scales both equal to PrW + ^, the cross sections for the three color-octet channels 
will be lowered down with the normalization factor changing from 1 to about 0.6 both for 
J/ip and ip'. But this change has little effects on the shapes of the curves in Figs. 1-2. So, 
our conclusion for the ^Sq dominance will not change with the scale changes. Another 
uncertainty may come from the parameterizations of the unintegrated gluon distributions, 
for which, however, it is difficult to give an estimate because we have little knowledge about 
it at present. On the other hand, we note that the parameterization provided by ||T^ are 
obtained from an excellent fit to F2{x,Q'^) in a very large window of x and Q^ in the kt 
factorization approach, which may be viewed as a reliable determination of these quantities. 
In conclusion, in this paper we have calculated the J/ip and ip' production rates and 
polarizations in hadroproduction at the Tevatron in the kt factorization approach. Various 
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color-octet contributions including ^5*0 , ^Pj , and ^S*} as well as contributions from Xcj 
decays are considered in this approach. We find that in a rather wide range of the trans- 
verse momenta of J/ip and ip', the production rates could be dominated by the ^Sq channel, 
and the predicted polarizations are in rough agreement with the preliminary Tevatron data. 
However, to draw the final conclusion on heavy quarkonium production mechanism in high 
energy hadron collisions, we still need more accurate experimental data on large pt pro- 
duction cross sections and polarizations. We also need more solid theoretical calculations 
such as NLO corrections to the color-octet processes discussed above, and a more reliable 
estimate for the so-called unintegrated gluon distribution function. It is worthwhile to make 
further investigations along this direction. 

While the calculation in this paper was finished we noticed a paper by Ph. Haegler et al. 
10| , who studied the J/ip production rate in the kt factorization approach with color octet 



channels. Their result is similar to ours for the calculation of cross sections. 
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Figure Captions 

FIG. 1. The pt distribution of direct J /ip production cross section at the Tevatron in the 

kt factorization approach. The sohd hne is for ^5*0 contribution, the dashed hne for ^Pj , 

and the dotted-dashed hne for ^S{ . 

FIG. 2. The same as Fig. 1 but for if)' . 

FIG. 3. Prompt if)' polarization at the Tevatron. 

FIG. 4. Prompt J /ip polarization at the Tevatron. The thinner lines are for the direct 

polarizations, and the thicker lines for the prompt polarizations. 
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